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GLOBAL WEAK SOLUTIONS FOR KOLMOGOROV-VICSEK TYPE 
EQUATIONS WITH ORIENTATIONAL INTERACTIONS 

IRENE M. GAMBA AND MOON-JIN KANG 


Abstract. We study the global existence and uniqueness of weak solutions to kinetic 
Kolmogorov-Vicsek models that can be considered a non-local non-linear Fokker-Planck 
type equation describing the dynamics of individuals with orientational interactions. This 
model is derived from the discrete Couzin-Vicsek algorithm as mean-field limit mm. 
which governs the interactions of stochastic agents moving with a velocity of constant 
magnitude, i.e. the the corresponding velocity space for these type of Kolmogorov-Vicsek 
models are the unit sphere. Our analysis for estimates and compactness properties 
take advantage of the orientational interaction property meaning that the velocity space 
is a compact manifold. 


1. Introduction 

Recently, a variety of mathematical models capturing the emergent phenomena of self- 
driven agents have received extensive attention. In particular, the discrete Couzin-Vicsek 
algorithm (CVA) has been proposed as a model describing the interactions of agents moving 
with velocity of constant magnitude, and with angles measured from a reference direction 
(see [IlIlIIBlEe]). 

This manuscript focuses on analytical issues for the kinetic (mesoscopic) description 
associated to the discrete Couzin-Vicsek algorithm with stochastic dynamics corresponding 
to Brownian motion on a sphere. More precisely, we consider the corresponding kinetic 
Kolmogorov-Vicsek model describing stochastic particles with orientational interaction. 


dtf + io- V^/ = -V^ • ifFo) + 

Fo{x, uj, t) = iy{uj ■ n(/))(/d — a; (g) uj)Q{f), 

n{J){x,t) = J{f){x,t)= [ K{\x - y\)u]f{y,u:,t)dydu], 

f{x,u},0) = fo{x,U!), X G U, CO £ t > 0, 

where / = f{x,co,t) is the one-particle distribution function at position x £ U, velocity 
direction w G and time t. The spatial domain U denotes either or T'^. The 

operators and denote the gradient and the Laplace-Beltrami operator on the sphere 
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gd -1 respectively, and // > 0 is a diffusion coefficient. The term Fo{x,uj,t) is the mean-field 
force that governs the orientational interaction of self-driven particles by aligning them with 
the direction S that depends on the flux J{x,t). 

This mean-field force is also proportional to the interaction frequency v. Its reciprocal 
represents the typical time-interval between two successive changes in the trajectory 
of the orientational swarm particle to accommodate the presence of other particles in the 
neighborhood. The function K is an isotropic observation kernel around each particle and 
it is assumed to be integrable in M. 

Following Degond and Motsch in m, the interaction frequency function u is taken to be 
a positive function of cos 9, where 9 is the angle between uj and fl. Such dependence of i' 
with respect to the angle 9 represents different turning transition rates at different angles. 
Hence, the constitutive form of such interaction frequency u{9) is inherent to species being 
modeled by orientational interactions. As in m, we assume that i'{9) is a smooth and 
bounded function of its argument. 

The kinetic Kolmogorov-Fokker-Planck type model with orientational interactions 
was formally derived in m as a mean-field limit of the discrete Couzin-Vicsek algorithm 
(CVA) with stochastic dynamics. There, the authors mainly focused on the model (II.ip 
with the following local momentum J instead of J 

(1.2) n{J){x,t) = J{f){x,t)= j ujf{x,u},t)du}, 

where J was derived from J in (II.ip by rescaling the kernel K in time and spatial vari¬ 
ables. Such scaling describes dynamics for solutions to (11.11) at large time and length scales 
compared with scales of the individuals. 

In the current manuscript, we focus on existence and uniqueness properties of solutions 
to both models, with J(/) as defined in (II.ip and with J{f) as defined in m- 

In fact, since J with the kernel K = 5o(Dhac mass) is exactly J, it is enough to show 
global existence and uniqueness of weak solutions to models (II.1|) in an appropriate space, 
to be specified in Section 2. These results are easily applied to J as in (11.21) . 

The classical Vicsek model have received extensive attention in the last few years con¬ 
cerning the rigor of mathematical studies of its mean-field land hydrodynamic limits as 
well as phase transition development. More specifically, Bolley, Cahizo and Carrillo have 
rigorously justified a mean-field limit in [2] when the force term acting on the particles 
is not normalized, i.e., vQ.{x,t) replaced by just J{x,t) in force term F),. This modifica¬ 
tion leads to the appearance of phase transitions from disordered states at low density to 
aligned (ordered) states at high densities. Such phase transition problem has been studied 
in pill El Emu [IB]. In addition, issues on hydrodynamic descriptions of classical Vicsek 
model have been discussed in PEiiiniiiiiiiaiisi. We also refer to laElllg] concerning 
related issues. 

Up to date, there are few results on existence theory of true kinetic descriptions. Frouvelle 
and Liu m have shown the well-posedness in the space-homogeneous case of dOD with 
the regular force field {Id — uj®uj)J instead of {Id — uj®uj)0,{J). There, they have provided 
the convergence rates towards equilibria by using the Onsager free energy functional and 
Lasalle’s invariance principle, and their results have been applied in [8]. Very recently. 
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Figalli, Morales and the second author m have shown the well-posedness in the space- 
homogeneous case of (jl. 2 j) . and the convergence of solutions towards steady states, based 
on the gradient flow approach (see for example nam]). 

On the other hand, the authors in [2] have shown existence of weak solutions for the space- 
inhomogeneous equation for a force field Fo given by the difference between spatial convolu¬ 
tions of mass and momentum with a bounded Lipschitz kernels K, namely uK*xp — K*xJ, 
instead of vQ. as considered in this manuscript. Such a choice of force field has a regularizing 
effect for spatial variable compared to our case which deals with stronger non-linearities. 

This manuscript is mainly devoted to showing the existence and uniqueness properties 
of weak solutions to the kinetic Kolmogorov-Vicsek type model (II.111 . A difficulty in our 
analysis arises from the fact that ff(J) in the alignment force term of (jl.lj) is undefined at 
J(/) becomes 0. So we restrict the problem of finding global weak solutions to (II.ip to a 
subclass of solutions with the non-zero local momentum, i.e. J{f) 7 ^ 0 . 

In the next section, we briefly present some known results for kinetic models with ori¬ 
entational interactions, 0 and which give a heuristic justification for the a priori 

non-zero assumption on J{f) to be stated in our main result. Section 3 presents a priori 
estimates and the compactness lemma, which play crucial roles in the main proof of exis¬ 
tence of weak solution in the next section. Section 4 deals with the construction of weak 
solutions to (diD by means of first, introducing an e-regularized problem, for an arbitrary 
parameter e > 0 modifying the alignment force n(J) uniformly bounded in e. We then 
solve the e-regularized problem of (jl.ll) constructing a sequence of functions {/n,e}n>i that 
converges to the solution /g. Finally we show that, in within the class of solutions satisfying 
J{f) > 0, there is a subsequence converging to /, solving (II.ip . Section 5 is devoted to 
the proof of the uniqueness of weak solutions in a periodic spatial domain U = under 
the additional constraint J{f) > a > 0. 

2. Preliminaries and Main results 

In this section, we briefly review how the kinetic Kolmogorov-Viscek equations, (II.ip and 
(jl.2j) can be formally derived from the discrete Couzin-Vicsek algorithm model [TO] with 
stochastic dynamics. Then we provide our main result and useful formulations. 

2.1. Kinetic Kolmogorov-Vicsek models. Following [lOj. the kinetic Kolmogorov-Vicsek 
model considered in (11.111 is derived from the classical discrete Vicsek formulation modeling 
Brownian motion of the sphere given by the following stochastic differential equations 
for 1 < z < V, 


( 2 . 1 ) 


dXi = ujidt, 

duji = {Id — uji® uji)v{uji ■ Qi)^idt + ■\/^{Id — uji® uji) o dBl, 

— I j I ) Ji — 'y ^ ■ 

' i, \Xj-Xi\<R 


Here, the neighborhood of the i-th particle is the ball centered at Xi G with radius 
R > 0. The velocity director w* G ^ of the z-th particle tends to be aligned with the 
director flj of the average velocity of the neighboring particles with noise Bl standing for 
N independent standard Brownian motions on with intensity v^2)u. Then, its projection 
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{Id — Ui® (jJi) o dBl represents the contribution of a Brownian motion on the sphere 
which should be understood in the Stratonovich sense. We refer to m for a detailed 
description on Brownian motions on Riemannian manifolds. We note that the first term 
in duji is the sum of smooth binary interactions with identical speeds, whereas there is 
no constraint on the velocity in the Cucker-Smale model [5]. In addition the interaction 
frequency (weight) function u{uji-lli) depends on the angle between Wj and Hi, parametrized 
by cos 6i = uji ■ Hi- 

From the individual-based model m, the corresponding kinetic mean-field limit (II.ip 
was proposed in mm, as the number of particles N tends to infinity. Notice that /i 
in dni) corresponds to the diffusive coefficient associated to the Brownian motion on the 
sphere 

The reduced model dni) with the modified definition of setting J = J as in ()1.2p was 
proposed in |10j by the following scaling argument. Considering the system dynamics at 
large times and length scales compared with those scales of individuals by the dimensionless 
rescaled varaibles x = ex, t = et with e <C 1 , it makes the interactions to become local and 
aligned the particle velocity into the direction of the local particle flux. This interaction 
term is balanced at leading order e by the diffusion term. 

Notice that H{f) in ( 11 . 11 ) is undefined when J{f) becomes 0. Because of this issue, we 
study in this manuscript the existence of weak solutions to for the subclass of solutions 
with non-zero local momentum, i.e. J{f) 7 ^ 0. As shown in [ID] , since oj is not a collisional 
invariant of operator Q, the momentum is not conserved. Thus, it is not straightforward to 
get J{f){x, t) 7 ^ 0 for all {x, t) from imposing non-zero initial momentum, i.e. J{f){x, 0) 7 ^ 0 
for all X. Moreover, there is no canonical entropy for the type of the kinetic equations as in 
(dU). Due to these analytical difficulties, we heuristically justify our constraint J(/) 7 ^ 0 
by observing equilibria of (jl. 2 p in the three dimensional case, which has been studied in 
m as follows. 

For the classification of equilibria in the d = 3 dimensional case, we recall the the Fisher- 
von Mises distribution, given by 


Mq{uj) 


/s2 exp do; 


exp 


cr(a; • D) 


for a given unit vector H E S^, where a denotes an antiderivative of u, i.e. ^(t) = 

Since v is positive, a is an increasing function and then Mq is maximal at w • D = 1 , that 
is for u! pointing in the direction of D. Therefore, H plays the same role as the averaged 
velocity in the classical Maxwellian equilibria for classical kinetic models of rarefied gas 
dynamics with velocities defined in all space. The diffusion constant /i corresponds to the 
temperature strength, which measures the spreading of the equilibrium state about the 
average direction H. The present model has a constant diffusion /r that is in contrast with 
the classical gas dynamics where the temperature is a thermodynamical variable whose 
evolution is determined by the energy balance equation. 

Using the Fisher-von Mises distribution, the operator Q and equilibria of (|1.2p are ex¬ 
pressed as follows. 


Lemma 2.1. [TO] (i) The operator Q{f) can be written as 





Qif) = 
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(ii) The equilibria, i.e. solutions f{uj) satisfying Q(f) = 0, form a three dimensional man¬ 
ifold £ given by 

£ = {pMn(w) I /O > 0, O e S^}, 

where p is the total mass and Q is the flux direetor of pM^{uj), that is, 


P = 


pMn{uj)duj, 


JjpMn) 

J{pMn) 


JipMn) := 


pM^{uj)ujdoj = pc{p)£l, 


with 

cos 0 exp ^ gin 8d6 

/J'exp sin0d0 

We note that c{p) —)■ 1 as /x —)> 0, and c{p) —)• 0 as // —>■ oo. This means that the local 
momentum J{pMq) of the equilibrium solution / = is not zero as long as the diffusion 
strength p is not sufficiently large compared to orientational interaction. Consequently, it is 
expected that moderate values of p would yield non-zero local momentum J(/) for solutions 
/ near the Von Mises equilibria. 


2.2. Main result. We state now the main results for global existence of weak solutions to 
equations dni). 

We hrst introduce the following notations for simplification. 

• Notation : We denote hy D -.= 11 x and by := Id — Lo iSi ui, as the mapping 

V !-)• {Id — uj ® uj)v is the projection of the vector v onto the normal plane to to. 

• Hypotheses {H) : As stated earlier, we assume that iy{-) is a smooth and bounded 
function of its argument and K{\ ■ |) G L^{U). Moreover, in order to avoid Il{f) to be 
undefined, we impose a priori assumptions stating that the weak solutions / of (11.11) belong 
to an admissible class 

(2.2) A:={f\ J{f){x, t)fl0, yxeU,t> 0}. 

Theorem 2.1 (Existence for spatial domains U, being either or T'^). Assume 
ifH) and fo satisfies 

(2.3) h£{L^nL^){D) and /o > 0. 

Then, for a given T > 0, the equation (HTD has a weak solution f, which satisfies 

/> 0 , 

/GC(0,r;L'(H))nL“(Hx (0,r)), 
v^/gl2(hx (0,r)). 


(2.4) 
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and the following weak formulation: for any (f> £ C^{D x [0,T)), 


[ [ fdt4> + fu: ■ Va;(/> + fFo ■ V^(/> - fiVujf ■ Vujcfdxdcods 
Jo Jd 


(2.5) 


+ / /o'(’(0, ■)dxdoj = 0, 

Jd 

Fo{x,uj,t) = u{uj • n{f))F^±n{f). 


Moreover, the weak solution f satisfies the following estimate 

(2-6 ) ||/||l°°(0,T;Lp(D)) + -^II^‘^/^IIl2(Dx(0,T)) - II/o||lp(D), 


for any 1 < p < oo, and 



Remark 2.1. The proof of Theorem \2.1\ is based on energy methods, where the diffusion 
term pA^^f plays a crucial role. Yet the strength p > 0 does not essentially affect the proof 
of existence. Therefore, without loss of generality, from now on we set p = 1. 


We next present uniqueness of weak solutions being constructed in Theorem l2.ll only for 
periodic domains U = T^, together with the following subclass 

:= {/ I 3 a > 0 s.t. |J(/)(x,t)| > a, Y{x,t) £ T'^ x (0,T)}, 

which is more restrictive than (j2.2p . Indeed this class corresponds to the subclass of weak 
solutions to the initial value problem (II.ip . with uniformly bounded below speed when 
solved in a spatial torus domain. 

Theorem 2.2 (Uniqueness for periodic spatial domains T^^). Assume {%) and (12.3p . 

Then for a given T > 0, the periodic boundary problem of o has a unique weak solution 
f in the subclass Aa- 


Remark 2.2. Our proof for uniqueness takes advantage of a uniformly positive lower bound 
a of J{f) in order to control YL{f), consequently restrict to periodic domain T'^. Indeed, 
imposing that J{f) > a > 0 for all x £ results in an infinite mass fdxdoj = oo, 

due to 


oo = 


[ J{f)dx < [ 

JR'* JR'^xS 

</ 


\K 


\K f\dxdu} 

*x f\\L^(Rd.-)duj 


/ fdxdw, 


2.3. Formulas for Calculus on sphere. We start recalling some useful formulas on the 
sphere which are extensively used in this paper. 

Let F be a vector-valued function and / be a scalar-valued function. The following formula, 
as a analogous of the integration by parts, holds 

(2.8) [ fV^-Fdw = - I F- (V^/ - 2uj f)du:. 

J§d-1 J§d-1 

By the dehnition of the projection operator it follows that 

F^±uj = 0, P^xV^/ = V^/, 

P^xU • V = P^jXU • u, 


(2.9) 
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for any scalar-valued function /, and vectors u and v. 


In addition, for any constant vector v G we have 


( 2 . 10 ) 


V^(cj • v) = 

• (Pa;-L^^) = -{d-l)uj ■ V. 


These formulas can be easily derived classical calculus on spherical coordinates (see mm)- 


3. A PRIORI ESTIMATES AND COMPACTNESS LEMMA 

The following Lemma provides a priori estimates in L°°(0, T; L^iU)), 1 < p < oo for solu¬ 
tions to the initial value problem for the kinetic equation below. The subsequent Lemma [3.2l 
provides a compactness tool needed for the existence result proof of Theorem 12.11 

Lemma 3.1. Assume that /o satisfies (12.3p . and f is a smooth solution to the equation 

^3 dtf + oo- VJ = -V^ • (/p(a; • L!)P^xL!) + A^f, 

f{x,uj,0) = fo{x,uj), 

where Q, : U x M+ is a bounded vector-valued function of {x,t). 

Then, for any 1 < p < oo, 

‘lin — 1 ) p - CT p 

(3-2) ||/||l,°°(0,T;I,P(D)) H- - - \\'^‘-}f^\\h{Dx{0,T)) - ^ ^"Ml/o IIlp(D) • 


In particular, ifp = oo, then 



Proof. First of all, for any 1 < p < oo, it follows from (II.1|) that 

f [ f^dxduj = —p I f^~^'Vu}-{fn{io-Q)F^±Il)dxdu}+p [ f^~^fdxdco 
dt Jd Jd Jd 


(3.4) 


= : h+h. 


Using formula (|2.8p and w • V^if = 0, we have 

I 2 = -p{p - 1) / • V^fdxdu + 2p [ ■ V^fdxdio 

JD Jd 

= f iv^f^'^dxduj. 

P Jd 

Next, by formula (I2.10|) . the term Ii from (13.4p is estimated as follows 
h = —p j ■ P^^xH + fvfuj ■ n)|P^xnp — {d — l)/p(a; • fl)w • P^dxduj 

< p\\n{u} ■ Q)\\l°° [ f^~^\Vu}f\dxdu}p\\n'{uj ■ [ f^dxduo 

Jd Jd 

-\-p{d — l)\\v{ijj ■ Pt)\\Loc, f f^dxduj. 

Jd 
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In addition, using Holder’s inequality, the first integral in the right hand side above can 
be estimated by 

J F~^\Vujf\dxdu} < (y Fdxdu?j ^ j F~^\Vi_,f\‘^dxdu?j 

= - ^y f^dxdio'j ^ (y l^ujf^fdxdu 


Then, we have 


h < 


2{p - 1) 


P 


|Va;/2 1'^dxdui + C{- 


P 


ID 


P 


I 


pdxduj. 


Finally, combining the estimates above for both Ii and I 2 , we get 
d r , 2{p-l) 


, [ Fdxdu + 
dt Jd 


P 


f I'^ujf^l'^dxduj < C (— - \-p) [ Fdxdco, 

Jd P-'^ Jd 


which yields a Gronwall type inequality 

^||/||lp(d) < 

Therefore, 

ll/llL°°(0,r;LP(D)) < 11/0II LP(_D), 

which implies the estimate in (|3.2p . Hence, taking p —>■ 00 , yields the L°° bound (j3.3p . □ 

Remark 3.1. The boundedness of the alignment vector TL is essential for the proof of Lemma 
\3.1l and the a priori estimates and (j3.3p still hold for Q = Ll{f) bounded for any f. 


(3.5) 


The following lemma provides the compactness property that ensures the strong con¬ 
vergence of solutions to the initial value problem associated to linear equation (j3.5ll . Such 
strong compactness property relies on the boundedness of both the force term and velocity 
space (notice that the velocity variable would be unbounded, we would have to use the cele¬ 
brated velocity averaging lemma [221125]). As mentioned earlier, the compactness property 
obtained from next lemma is crucial for the existence proof of Theorem 12.11 

Lemma 3.2. Assume that /o satisfies (j2.3p . and fn is a smooth solution to 

dtfn + W • V^fn = -V^ • (^fnl^iuj ■ F„)P^±F„^ A^/„, 

fn{x,U},0) = fo{x,Uj), 
where Fn : U x M+ is a given function of {t, x). 

If the sequence (Fn) is bounded in L°°{U x (0, T)), then there exists a limit function f such 
that, up to a subsequence, 

fn^f as n ^ 00 in L^{D X {Q,T)) L?{U X , 1 < p < 00 . 

Moreover, the associated sequence 

Jn-.= / K{\x - y\)ujfn{y,oj,f)dydu) 

Jd 

strongly converges to the corresponding limit J in IF{U x (0, T)) where 


J ■■= j K{\x-y\)ujf{y,uj,t)dyduj. 


Id 
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Proof. Since the sequence (F„) is bounded in L°^{U x (0, T)), there exists F e L°°{U x (0, T)) 
such that, up to a subsequence, 


(3.6) 


F weakly — * in L°°{U x (0,T)). 


Let / be a solution of (I3.5p corresponding to the limiting F. Then, the following identity 
holds 


dtUn - /) + W • Vx{fn -f) = -V^ • ((/n - fHcO ■ 


(3.7) 


- V, 


- V,. 


(/(Ku; • TO - Z^(a; • F))P,xF0) 

■ T)P,x(T„ - T)) + AMn - /)• 


Next, for any fixed p G [l,oo), multiplying the above equation by p{fn — fT ^ and 
integrating over D yields the identity 


fL 

dt 


'D 


ifn - ffdxduj 


(3.8) 


+ 


pJo^fn- • [ifn - • iO)P^xio)dxda; 

pJo^fn- • [fiHu; • Fn) - v{uj • T))P^xdxdo; 

pj^ifn- • [MuJ • T)P^x(F„ - F)yxduj 

P [ (fn- f)^~^A^{fn - f)dx(hj 
JD 


—: ffi + 772 + iTs + TTj.. 

We first estimate the term J'l using the same arguments as the ones used in Lemma l3.II 
in order to estimate Ji. Indeed, 

Jl = -pjyfn - • Fn)^Mn " /) ' P..xiO + {fn - /)T(u; • TOlP^xi^P 

- {d- l){fn - f)v{!^ • Fn)uj • Fr^dxdw 
< p|0(w • iO)||i,oo||F„||Loo [ {fn - fT~^\'^uj{fn - f)\dxduj 

JD 

+ p\\u'{uj ■ Fn)\\L°°\\Fn\\loo [ {fn - ffdxduj 

JD 

+ p{d - l)\\l'{uj ■ Fn)\\L<^\\Fn\\L^ [ {fn - ffdxcLj. 

JD 

<— - [ \Vaj{fn- f)^fdxdCL! + [ {fn- fYdxduj. 

P Jd P-^ Jd 

Similarly, 774 is also estimated as done for I 2 in the proof of Lemma 13.11 

4{p - 1 ) 


774 = - 


P 


[ \^Mn-f) 

JD 


— f)2 rdxdu}. 
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Hence, gathering these two last estimates, identity (|3.8p yields the following estimate 

4 [ ifn- fTdxdu <C [ {fn- fYdxduj —— [ \Vu,{fn-D^l'^dxduj + J 2 +Js- 

dt Jd Jd P Jd 

Next, since /„ = / at t = 0, then applying the Gronwall’s inequality to the above 
inequality, it holds that for any 0 < t < T, 


[ {fn - ffdxdoj + 

Jd 


2{p - 1 ) 


P 


0 JD 


|Va;(/n - /)2 pdxdwds < [ {J 2 + J^){s)ds. 

JO 


The terms J 2 and can be rewritten using the calculaus on the sphere formulas (I2.10p 
as follows. First, note that the term J 2 satisfies the identity 

J2 = -P [ {fn - ff-^ \{H^ • Fn) - V{UJ • F))V^/ • 

Jd l 

+ f{u'{u: • Fn)Fn - F{u • F)F) • P^xF„ 

- {d — l)f{u{uj ■ Fn) — ^{u) ■ F))oj ■ Fn dxdio 

= -p [ {fn - /f W{UJ • F:)u, ■ {Fn - F)V^f • P^xF, 

+ f(u'{uj • Fn){Fn -F) + v"{u: • Fr)a; • (F„ - F)f) • P^xF„ 

- (d - l)fv'{uj ■ Ff)u ■ {Fn - F)uj • Fn] dxduj 

= -p [ {fn - f)^-^ W{uj • F:)V^f • P^xF^o; 

+ fu'iio • F„)P^xF„ + • F^F • P^xF^a; 

- (d - l)fF{uj • F^)a; • F^wj • (F„ - F)dxdu}, 

where F* and F** are some bounded functions due to the mean value theorem property, 
depending solely on the known bounded functions Fn{x,t) and its limit F defined in (13.6p . 
Similarly, also by the identities in (|2.10l) . the term jTa satisfies the identity 

J3 = -P [ {fn - /f ko; • F)V,/ • P,x(F, - F) + fu'{co ■ F)P,xF • P,x(F„ - F) 

— (d — l)fi'{ui ■ F)ui ■ {Fn — F) dxduj 

= -P [ {fn - f)^-^ UuJ • F)V^/ + fu'{u • F)P^xF - (d - l)fi^{u • F)a;l • (F„ - F)dxdu. 

Jd l -I 

Thus, we get the weighted estimate 

Wfn - /IIlp(z)) + ^ I'^Mn - /)2pdxda;ds 


< e‘ 


CT 


/ / $(x,r(;,s) ■ {Fn — F)dxdujds, 

Jo Jd 


(3.9) 
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where the weight function, given by 

ci>(x, w, s) = -p{fn - p(a; • F*)V^f • + fv'(u: • F„)P^xK + fv"{p: • • P^xF„a; 

- (h - l)/z.'(u; • F„*)u; • + v{u: • F)V^/ + fu'{u; • F)P^xF - (d - l)/i.(u; • F)uj . 

is shown to satisfy $ G L^{D x (0,T)). 

In order to show this assertion, first we show the uniform control property of both /„ and 
/ and their gradients. Indeed, by the uniform boundedness of (F„), applying the same 
estimates as in Lemma l3.II for both g = fn and /, respectively, we obtain 

ll5llL°°(0,r;LP(D)) < C'll/o||LP(i:))! 1 < p < oo, 

\Nujg^\\mDxiO,T)) < C'll/olliip), 1 <p < oo, 

where the positive constant C only depends on p and T. 


Next, by Holder’s inequality follows that 

rT 


lo Jd 


{fn - ff ^v^fdxdujds < ( y ifn- f)^dxdu?j ^ j {fn- fY ‘^\v^f\^dxduj 

C( y^(/^ + F)dxdw) ^ I l^dxda; 


1/2 


< 

<6-0, 


and 


J J {f^~ fy ^fdxdujds < ( y {fn - fYdxdu?j 

{Fn + 


P-1 

P 


pdxduj 


D 


< C 

<Co, 


D 


P-1 

P 


where the positive constant Cq depends only on ||/o||LP(ri)- 

Therefore, the weight function can be estimated by 

II‘1||l1(-Dx( 0,T)) < C'*(ll(/n — fY ^'^u)f\\F{Dx(0,T)) + ll(/n “ fY ^ f W F (D x (0 ,T))) 

< C^Cq, 

where the positive constant C* is given by 

C* = pd (^\\u'{u] ■ F*)||j;,oo + \\u'{lO ■ Fn)\\L’=°) + ' -^n*)l|L°° ||P||l°°^ ||Fn||L°° 

+ ||i^(a; • F)||loo + \\v'{uj ■ F)||loo||F||loo , 

which does not depend on n thanks to the uniform boundedness of the sequence F„. 

Hence, applying (13.611 to (13.91) . it follows that 

fn^f inFP(Fx (0,r)), 

V^/n^V^/ inL2(Fx (0,r)). 


(3.10) 
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Finally, in order to complete the proof of Lemma 13.21 it remains to show that (I3.10p 
implies the strong convergence of the associated sequence {Jn) = {J{fn)) towards J(/). 
Indeed, Minkowski inequality. Holder’s inequality and Young’s inequality yield 

rT 


\Jn — J\\lp{Ux{0,T)) — 


/I/ 


K *x ifn - f)u}diO 


0 JU 


(3.11) 


< 


f ([ [ \K*x{fn- 

^ Jo Ju 


-^Uo 

< C\\K\\L^(S.d)\\fn - /lll,P(Dx(0,T))) 


dxdsj^ 
^dio 


which completes the proof. 


□ 


4. Proof of Existence - Theorem 12.11 

The proof of Theorem 12.11 entices the construction of an iteration scheme that generates 
a sequence (/n), where /„ is a solution to the linear equation (j3.5p at n-th step, with 
Fn := evaluated at the (n — l)-th solution fn-i obtained in the previous (n — l)-th 

step. 

This first intuitive approach confronts a difficulty since such n-iteration scheme generating 
the sequence /„ does not secure the non-zero momentum \ J{fn)\ > 0, even if | J(/„_i)| > 0. 
In fact, if that would be the case, the term Q{fn) would be undefined and therefore we could 
not secure it is bounded. In particular, since the compactness properties of Lemma l3.II and 
Lemma 13.21 require a bounded force term ( in (13.11) and ()3.5p respectively) then, with with 
at least the available tools developed in this manuscript, it would not be possible to secure 
an existence of a solution fn+i for the next n-iterative step. 

Hence, a way to avoid this difficulty can be acomplished by the use of an e-regularization 
approach by adding an arbitrary e > 0 parameter to the denominator of fl(/„), for all 
n G N. Such regularization generates a double parameter (e, n) sequence of solutions fe,n 
that it is shown to satisfy the property | J(/e,n)| > 0 for all n G N, uniformly in e > 0. 


4.1. The e-regularized equation. The e-regularization approach consists in solving the 
non-linear problem (jl.ip by adding e > 0 to the denominator of 11(/), 


(4.1) 


dtfe + u; • Vxfe = -V^ • 




Je{fe){x,t) = 


\Je{fe){x,t)\ +e' 

fe{x,UJ,0) = fo{x,Uj), X G [/, W G t > 0. 


/GxS'*-! 


K{\x 


y\)^^fe{y,(^,t)dyduj 


This new non-linear e-problem is then solved by generating a sequence of solutions fs^n 
to (13.5p with a bounded ■= f^e(/£,n-i) for the previous iterated solution fe,n-i- 
In the sequel, we show first that is possible to construct a sequence of solutions 
converging to in LP{D X (0, T)) (H L‘^{U X (0, T); FJ^, 1 < p < oo for any e > 0, so 
the results remains true in the e —>■ 0 limit. 

The details of this procedure are as follows. 
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4.2. Construction of approximate solutions. The construction of an (e, n)-sequence of 
approximate solutions to the non-linear e-regularized equation (14.Ij) is now done by the 
following iteration scheme. For any fixed e > 0, set := fo{x,u}) to be the initial 

state associated to dni). Then, define as the solution of the following linear initial value 
problem 


dtfe,l + W • Vxfe,l = -Vuj • 

Je,o{x,t) = j K{\x - y\)ufefl{y,u},t)dydu} 

fs,l{x,UJ,0) = fo{x,uj). 

Inductively, each fe^n+i is define to be the solution of the following linear initial value 
problem 

dtfe,n+l T ^ ^xfe,n+l ^u) ' ^u)fe,n+l: 

(‘^•2) Og^(x,f)= , Je,nix,t) = [ K{\x - y\)u}fe,n{y,‘x,t)dyduj 

fe,n+l{x,UJ,0) = fo{x,Uj). 

The justification for unique solvability of the (e, n)-approximate initial value problem 
for n > 1, follows form the next lemma. 


Lemma 4.1. For any T > 0, e > 0, n > 1, assume that is a given integrable function 
and /o satisfies (m. Then, there exists a unique solution fe^n+i > 0 to the equation (021) 
satisfying the L'^-estimates: for any 1 <p < oo, 


(4.3) II/, ,n+l\\L°°(0,T;LP{D)) + 


2(P - 1) 


P 




£ 

2 

£,n+l 


1 1 ^ CT^ 

lL2(nx(0,T)) - ^ 


LP{D), 


and 


(4.4) ||/£,n-n||L°°(Dx(0,T)) < e'^^||/o||L°°(D). 

The proof of Lemma |4.II follows the same argument as Degond’s proof in [^. We include 
its proof in Appendix for the reader’s convenience. 


4.3. Passing to the limit as n —)■ oo. The convergence of fe^n towards some limit function 
/,, which solves the regularized equation ()4.1h . is secured by the following proposition. 

Proposition 4.1. For a given T > 0 and arbitrary e > 0, if fo satisfies ()2.3I) . then there 
exists a weak solution fe>0to equation 0T]) satisfying the -estimates: for 1 <p < oo, 

(4.5) ||/£||l°°(0,T;Lp(D)) + -^l|Va;// IIl2(£,x(0,T)) - II/o||lp(D), 


and 

(4.6) 


II/£||l°°(Dx(o,t)) < e'^'^||/o||Loo(£)). 
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Proof. Since the sequence defined in (14.2p is bounded in L°°{U x (0,T)), we use 

Lemma 13.21 with = ^s,n- Thus, there exists a limit function such that, up to a 
subsequence, 

fe,n^fe as n ^ oo in X (0,r))nL2([/x (0,T);Fi(S'^-i)), 






Js as n —)• oo in L^{U x (0,T)), 


that yields 


Indeed, 






Tel + e 


as n —)■ oo in L°°(0, T; L^(L))). 


/{/ 


irie.n - 


lU 


^{Je,n — Je) + |Te|(Te,n ~ Tg) + Je(| Je| — |Je,r 


(|Te,n| + £)(|Te| + s) 


dx 


< 


^iJe,n Te) + I Je| (Je^n Te) + Je(|Je| iTe^ril) 


|Te| + e 


eP Ju 

< C{e) {\Je,n - J^r + \Je,n - Je\^ + ||Ts,n| " \Je\\^)dx 

< C{S) [ I Je,n - Je\Pdx. 

■JU 

Therefore, the limit /e satishes the following weak formulation of (|4.ip : for all (p G 
C^iDx[0,T)), 

ft f 

fedtp + feUJ ■ Vx(p + feFe ' - V^ife ' Vuifpdxdujds 



0 JD 


+ 


f fo4>{0, ■)dxdu = 0, 
JD 

JeiFi 


Fe = z/(a; • t) = 


|Je(x,t)| +e' 

In addition, using Lemma [3.11 together with the boundedness of fig above, the estimates 
from (14.51) and (|4.6I) follow. □ 


4.4. Passing to the limit as e —^ 0. The proof of Theorem 12.II is completed after showing 
the convergence from (14.11) to dni) as 0 < e —)• 0, in the weak sense. In fact, it is enough to 
show such limit for any convergent sequence 0 < —>■ 0. 

First, consider a sequence 


Fk := 


T£fc 

Jck I + 


T£fc 


ujfe^du. 


Since such sequence is bounded in L°°{U x (0,T)) uniformly in £k, Lemma 13.21 can be 
applied, so there exists a limit function / such that, up to a subsequence. 


/e, ^ / as oo in LP{D x {{),T)) x (0, T); 

Jgj, —)• J as A: —)• oo in ]J’{U x (0,r)). 


(4.7) 
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Next, in order to see that / is the weak solution to (II.ip it is enough to show that / 
satishes the weak formulation ()2.5p as a limit of the following formulation for (j4.ip : 

rt 


fl 

Jo Jd 


fekdt4> + 


OJ ■ Jf 




Jet \ + 




Efc 


Jet \ + ^k 


- V^/e, • V^cpdxduds 


+ 


[ /o0(O, ■)dxdu} = 0, 
Jd 


for any cf) e C^{D x [0,T)). 

By the convergence of in (14.71) . clearly all linear terms in the above formulation 
converge to their corresponding terms in (j2.5p . On the other hand, the convergence of the 
nonlinear term requires further justification provided in the following Lemma. 


Lemma 4.2. Assume \ J{x,t)\ > 0 as in 




U! ■ J, 




Jeu \ + ^k 


P,. 


J. 


£k 


\Jek \ + 


. Then, as k 
V i^(j)dxdujds 


ff 

Jo Jd 




oo, 


OJ ■ J 


. —— • Vujfpdxdojds. 

\d\ 


Proof. By the properties (12.9p of calculus on the sphere applied to the projection opera¬ 
tor P^x is the identity operator acting on gradient functions of the sphere that is 

P^^x • holds for any test function of tc G Then, the limit as k ^ oo 

in (14.8p is identical to show the analog limit for the formulation without the projection 
operator. That is, for k ^ oo 
(4.9) 

W Jet ^ Jet V7 . f f 

\J\ J |J| 




-'Efc 


deA ^k' I I T 


V i^(j)dxdojds 




We hrst control the integrand in (j4.9p using the estimates (14.6p and boundedness of v, 
so that there is a uniform constant C such that 


fsk'^ 


OJ ■ Jf 


£k 


J. 


Sk 


Jetl J ^k' \ J£k\ + ^k 


L°°(Dx(0,r)) 


P \\f£k\\L^{Dx{0,T))\\k'\\L'=° < C, 


which implies, for some F, that 


fet’^ 


OJ ■ Jf 




J. 




deA ~k ^k^ \deA ~k ^k 


F weakly — * in L°°{D x (0, T)). 


Then, it remains to show that 
F = •' 


|J| J \ J\ ’ ^ xU xS<^ 1 I |J(x,t)| > 0}. 

In order to obtain this last identity, we consider the bounded set 

:= {{t,x,oj) G (0,r] X Br{0) x | |J(x,t)| > J}, 

where R and 6 are any positive constants, and Br{0) denote the ball in U, with radius R, 
centered at 0. 


Since f and —)• J a.e. on Xr^s by (|4.7p . then by Egorov’s theorem, for any 

r] > 0, there exists a C Xr^s such that < i] and 


fek ^ di 


Sk 


J inL'^PP,), 
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and so, for sufficiently large k, 


\Jskix,t)\>- foT{x,t)eY^. 
Therefore, the L^(Y^) e-convergence follows from 



J 



\Jek 1 + ^k 

1^1 


(Mefc + ^k)\J\ 




^2 |J|(J,^-J) + J(|J|-|J,J)-e,J 


1^1 


L^iYv) 




that yiels 


fsk’^ 


OJ ■ Jt 


Efc 






^€k\ 






+ 


fek 


Ul ■ Jf 


£fc 


Jek \ + Efc 
OJ ■ J 


— V 


J. 


£fc 


OJ ■ J 




L^iYv) 


J. 


£k 


J 


l*^! ^ ^ \Y£k |<^| 
^ Cll/efc ||l°°(||z^^||l°° + l|2^l|i°°) 


L-(y,) 


+ 


{h,-M 


J. 


£k 


'^si. I T ^ \'^£k\ Y £k 




•^£k 

J 

\Jek 1 + ^k 

\J\ 


L“(n,) 


+ C'll/efc - /llL°°(Y„)lkl|L°° -> 0 . 


Hence, the following identity holds 

and, since 77 , ii and <5 are arbitrary, taking 77 , J —)• 0 and ii —)• 00 , it follows 

e (0,T] x [/ x S'^"^ | |J(x,t)| > 0} . 

which completes the proof of Lemma 14.21 □ 

Finally, thanks to lemma 021 and ()4.7I) . it follows that / satisfies the weak formulation 
(12.51) . In addition, estimates ()2.6p and (j2.7p follow directly from p4.5p and (14.6p . respectively. 
Therefore, the proof of Theorem 12.11 is now completed. 


5. Proof of Uniqueness - Theorem 12.21 

The uniqueness argument is considered in the subclass Aa of weak solutions constructed 
in Theorem 12. II Let / and g be any weak solutions to the initial value problem (11.11) in Aa- 
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A straightforward computation yields that 

^4/ \f - gl^dxduj + [ \V^{f - g)\'^dxdu} 


(5.10) 


[ if - g)^uj ■ (if - g)f{ud ■ n{f))¥^±n{f))dxdoj 
[ if - g)'^uj ■ • fl(/)) - u{uj ■ n{g)))F^±n{f))dxduj 

[ if - gWuj ■ (gj^i(^ ■ Ft{g))¥^± (fl(/) - Flig))) dxdoj 


— : Ji + J 2 + J 3 . 

Using the same estimates applied to J7i in the proof of Lemma l3.21 we can also estimate 

Ji<]-[ \V^{f - g)\'^dxcLj + C f \f - gl'^dxdu. 

4JT‘^xS‘*-i JTd^§d-i 

Next, J 2 and J 3 can also be estimated same approach from (13.111) in Lemma [321 to get 

ll'^(/) “ dig)\\L^(Td^ ^ “ 5'llL2(T'ixS‘*-i)' 

Moreover, since | J(/)| > a in the set Mq,, then the difference of alignment forces for any 
two weak solutions is controlled by 

\J{g)\{J{f) - J{g)) - J{g){\J{f)\ -\J{g)\) 

\nif)-n{g)\< - 


a\Jig)\ 


<-\J{f)-J{g)\ 


a 


that yields 

ll^(/) “ ^(fl')llL2(T‘*) — C'll/ — 5|lL2(']rdxS‘*-i)- 

Therefore, by property (12.4p for any weak solution, the control of term J 2 in (I5.10p follows 
from 

J2 = / V^(/ - g) • ¥^±Q{f)g{u{uj ■ Q{f)) - v{oj ■ Ft{g)))dxduj 


< lbllL°° ||i^' 

si/ 

4 Jfdxi 


’'IIloc / \S/ ^if - g)\\nif) - n{g)\dxduj 

i^if — g)\‘^dxdoj + C f \f — gl'^dxduj. 

1 JTdxSd-l 


Likewise, the control of the last term J 3 in (|5.10l) follows, since 

J 3 < ^ / \Vujif-g)\‘^dxduj + C [ \f-g\^dxduj. 

^ JTdy^gd-l JTdxSd-l 

Hence, gathering the above estimates and using Gronwall’s inequality, we have 
/ \f - gl'^dxdu < f \fo - go\‘^dxduj, 

Jjd^^d-l JTpdxSd-l 

which implies the uniqueness of weak solutions to the initial value problem pi.ip in Ao 
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6. Conclusion 

We have shown the existence of global weak solutions to problem (as well for J 

dehned as in (ll.2p l in a subclass of solutions with the non-zero local momentum. These 
solutions are unique on the subclass of solutions in the d dimensional torus whose mean 
speed is uniformly bounded below by a strictly positive constant. 

An important future work would be to remove our assumption on the non-zero local 
momentum. The main difficulty is due to the lack of momentum conservation for solutions 
to problem (II.ip . and canonical entropy associated to the equation in (II.ip . Thus, at this 
point, we have neither suitable functional spaces nor distances to study the behavior of 
solutions whose momentum may vanish locally. This difficulty is related to the issue on 
stability of solutions to m- Another future work is to extend the uniqueness result to the 
whole spatial domain 


Appendix A. Proof of Lemma 14.11 

For the notational simplicity, we omit the subindex n -|- 1 in (14.21) . Our goal is to prove 
existence of solutions / to the linear equation 

dtf + uj- VJ = • (/u(a; • + A^/, 

(^• 1 ) ^ = TTr%r^’ K{\x - y\)ujg{y,uj,t)dyduj, 

\J[x,t)\+£ Jd 

/(x,w,0) = fo{x,uj), 

where g is just a given integrable function. 

We begin by rewriting (jA.ll) as 

dtf + oo- VJ + v{u: • • V^/ 

(A.2) + fu'{uj • fi)|P,^xOp -{d- l)fu{uj ■n)uj-n- A^f = 0, 

f{x,uj,0) = fo{x,uj), 

where formulas (j 2 . 10 p on projections and calculus on the sphere were used. 

Next, taking f{x,u!,t) := f{x,uj,t) for a given A > 0, It leads to the modified initial 
value problem 

dtf + UJ ■ Va;/-)-i/’i • Vujf + 1/^2 Pi/’s)/- A^/= 0 , 

f{x,uj,0) = fo{x,u), 

where the functions 'ipi, i /’2 and 1/^3 are given by 

il!i{x,io,t) = v{uj ■ fi)p^xfi, 

'll^ 2 {x,UJ,t) = v'{uj ■ fi)|P^xfip, 
l/^ 3 (x, UJ,t) = —{d— \)v{uj ■ fijcj • , 

respectively. Now, since |fi| < 1 and the smooth function v is bounded, then ipi, 'tp 2 and 
"03 are also bounded. Therefore, by J. L. Lions’ existence theorem in [23], the existence of 









19 


a solution for (IA.3P follows from the same argument given by Degond in [6]. That means, 
equation (IA.3(i has a solution / in the space 


y := {/ G L^([0, T] X u- \ Otf + UJ ■ V,/ G L^([0, T] x U' 

Furthermore, by the Green’s formula used in [6], then the following identity holds, for any 

/ey, 

(A.4) {dtf + 0 J- Va^fJ) = i ^(|/(x,w,T)p - \f{x,u},0)\‘^)dxdoj, 

where (•,•) denotes the pairing of L^([0,r] x U] and L^([0,r] x U] . 

This identity (IA.4P is needed below to show uniqueness of solutions / in y as follows. 

Let / G y be a solution to (|A.3I) with initial data /o = 0. Then, by ()A.4p . it follows 


(A.5) 


0 = {dtf + UJ ■ Vccf + V’l • Vo;/ + {X + i>2 + V’a)/ - A^/, /) 


'D 


\f{x,uj,T)\^dxduj — - I Vuj ■ ipi\ffdxdu 


ID 


+ J {\ + ip2 + 'ilJ3)\f\‘^dxdu;+ I \Vojf\‘^dxdui 


Id 


^ ■ V’i||l°°([o,t]xD) — IIV’2||l°°([o,t]xD) — IIV’3||i,°°([o,r]xZ))) j I/I dxduo. 


Next, since 


= v'{uj ■ ii)|P^xOp — {d— ■ 0)a; • Cl, 

the term • V’l is bounded. Thus, choosing A such that 


(A.6) 


1 


A > -||Vtj • 0 i||l°°([O,T]xD) + ||02||l°°([O,T]xD) + ||0’3||l°°([O,T]xD)! 


then, estimate (|A.5p yields / = 0, which proves the uniqueness of the linear equation (IA.3p . 
Therefore, (IA.3h has a unique solution / G L^([0,r] x [/; 

Furthermore, since fo>0 and /o G L°°{D), by a similar argument as in (IA.5h . 

/>0 and /G L'^([0,r] X D). 

Indeed, using the following identity from [6] on any / G T, with /_ := max(—/, 0), 

{dtf + u} = ^^(|/-(x,w,0)|^ - \f-{x,u!,T)\^)dxdu}. 
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Then, since /_(a;,a;,0) = 0 when /q > 0, it follows 

0 = (<9t/ + w • Va;/ + '0i •V^/+(A + V'2 + V'3)/-A^/,/_) 

=-]: [ \f-ix,uj,T)\'^dxdu}+ ]- [ ■ 'ijji\f-\‘^dxduj 

^ Jd ^ Jd 

- / (A +-02 + - [ \V^f-\‘^dxdu; 

J D JD 

< — — • '0i||l°°([o,t]xD) — \\'4’2\\l°^{[o,t]xD) “ IIV'3||l°°([o,t] xD)) / If-l'^dxdoj. 


Using the same A as in (IA. 6 I) . yields /_ = 0, which proves / > 0. 

The same argument also deduces that 

II/IIl°°([0,T]xD) < II/o||l°°(D)- 

Finally, using the transformation f{x,uj,t) = f{x,io,t), the results hold for solutions 
of (jA.2h as well. In addition, since the / properties are invariant under such transformation, 
then the proof of existence is completed, and estimates (14.31) and (14.41) follow directly from 
Lemma l3.II together with boundedness of 0. 
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